Initial states of qubit— boson models leading to conserved quantities 
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There are initial conditions of composite qubit-environment systems which provide conservation 
of certain qubit observables but are not associated with a conservation law related to a symmetry 
of the system. This states are shown to form a subspace of the qubit-environment space of states. 
General construction of such states is presented and illustrated by two examples: The first one is 
exactly solvable Jaynes-Cummings model and the second is the multi-photon Rabi model. 



CO 

o 

CN 
C 

a 
m 

CN 



43 : 
a* 
^— > ' 



o 

CO 



X 



PACS numbers: 03.65.Yz, 03.67.-a 

Keywords: Open quantum systems, qubit-boson model, conservation law 



Introduction. Interaction of two quantum systems al- 
most always modify a set of quantum numbers suit- 
able for the description of non-interacting components. 
Generically, only in the absence of interactions the set 
of good quantum numbers of the composite system con- 
sists of good quantum numbers of its constituents. Let 
us consider the energy as an example. It is known that 
(Hamiltonian— type) interaction of two quantum systems 
causes an energy flow from one system to the other 
and/or vice-versa. The total system remains conserva- 
tive but the character of the subsystems changes qual- 
itatively. They become open [1]. Interestingly, there is 
a very specific situation, commonly refereed to as de- 
phasing (or pure dccohcrcnce) [2, 3], when the energy ex- 
change is absent during the entire evolution, regardless of 
the initial conditions. However, for such phenomenon to 
occur one requires a rather uncommon feature of the in- 
teraction resulting in the existence of conservation law(s). 
It is also the case if, instead of the energy, one considers 
any observable acting on the constituent of the total sys- 
tem. One cannot expect its conservation unless certain 
symmetry is present in the system. 

At this point one can pose a natural question: Is it pos- 
sible to design quantum evolution that guarantees con- 
servation of, at least, some observables regardless of the 
existence of obvious symmetries in the system? In this 
Brief Report we provide the answer in the case of two- 
level quantum systems interacting with the environment. 

The layout of this work is as follows. We begin by 
showing how a proper choice of an initial preparation 
of the composite qubit-boson system can assures no en- 
ergy exchange between the subsystems. In other words, 
we identify dcphasing-like behaviour in non-dephasing 
model. Next, in section II a general method is presented. 
Section III followed by the conclusions IV serves as our 
second example. 

First example. Let us start with a simple exactly solv- 
able Jaynes-Cummings model [4]: 

H = ljo z + va'a + (g*cr + ® a + go-- ® a') , (1) 



where a and are the creation and annihilation opera- 
tors of the bosonic field (environment) i.e., [a, a'] = Ib 
and o~ z denotes the one of the Pauli spin matrices with 
eigenstates |±), whereas v and w describe the qubit and 
field energies, respectively. The coupling constant g re- 
flects the strength of the interaction between the systems, 
while c±|=f) = |±) and <r±|±) = are ladder operators. 

Although the model (1) essentially originates from 
quantum optics [4], it has been studied in wide range of 
other branches of physics (see e.g. [5]), for almost half- 
century, within the broad variety of contexts [6]. There- 
fore, it serves as a good starting point for our considera- 
tions. 

In this example we are interesting in finding qubit- 
boson initial state(s) g, such that the qubit energy 
Eq = wTr (<7 Z ® lBg(i)) ~ (c z ) is conserved, or Eq = 
LoTr (a z p(t)) = est. where the reduced qubit dynamics 
p(t) = Tr B (e- fflt ge mt ) is used [7]. 

It is well known that for the Jaynes-Cummings model 
the total number of excitation N = a^a + a z is conserved 
which ultimately leads to the exact diagonalization of the 
model and in particular results in exact solvability of its 
dynamics. At this stage we want to emphasize that none 
of the above features is essential in the following analysis, 
yet involving N one simplifies calculations. 

We start by splitting the Hamiltonian (1) into two com- 
muting parts, H = Hq + V, where Hq = vN and 



V = Sa z + (g*cr + <S> a + go- ® a') 



(2) 



where 5 = u> — v is the detuning frequency. As the result 
the time evolution of the systems can be factorized such 
that 



U t = exp(-iH t) exp(-iVt) = U?V 



(3) 



Let us observe that in this example p(t) mimics the pure 
dephasing evolution if 



p(t) = a\+)(+\+c(ty\-)(+\ 

+ c(t)\+)(-\ + (l- a )\-)(-\, 



(4) 
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where a is a real constant and c(t) denotes a function of 
time such that \c(t)\ < 1. 
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At this point two questions can be addressed: Which 
initial states of the composite system, g, do guarantee the 
reduced dynamics (4)? Are they separable or entangled? 

In order to answer this questions we offer an explicit 
construction of such states. At face value subsequent 
steps may seem bizarre and to some extent artificial. 
Notwithstanding this, they can be performed, as it is 
discussed below, for a broad class of qubit-boson mod- 
els. As a first step we define g = |^ r )($|, where 

|*) = CV(|+)®|V> + |-)®X|^)), (5) 

with \tf>) an arbitrary state of the bosonic field. is a 
normalization constant and (^l^) = 1. For the sake of 
simplicity we absorb it into \tp) so that \\ip\\ = \C^\. 

Possibility of constructing dephasing states relies on 
the ability of finding a linear operator X satisfying 

g*XaX + 2SX - go} = 0, (6) 

or X(<5 + g*aX) = —5X + ga^ . As it will be seen, this 
choice of X allows for a specific time evolution of the 
above state, namely 

|*) -> |tf t ) = U° (|+) ® |^) + |-) ® X|^» . (7) 

According to (7), the dynamics of the entire system con- 
sists of two parts. The one governed by U° = exp(— iH t) 
is the free evolution. The other one, encoded in \ipt), is 
the remaining part resulting from the interaction and it 
needs to be determined along with X. This results in the 
following reduced dynamics 

<*(*) = &t\i> t ), 

c(t) = e- i2ut {ijjt\e iva]at Xe- iv ^ at \ijjt), 

which resembles the dephasing dynamics (4) if and only 
if the vector \ijjt) itself evolves unitarily. For this to 
hold true there must exist a Kamiltonian K such that 
\tp t ) = e~ iKt \ili). In order to see this is indeed the case 
note that according to (6) and (2), for K = 5 + g*aX one 
has 

V|tt) = |+) ® K|V> + |-) ® (-<5X + .gat) |V>) 
= |+)®K|V) + |-)®XK|V>), 

for every \tp). Thus, by replacing with K|t/>) one finds 
V 2 |*) = |+) ® K 2 |^) + |-> ® XK 2 |^) (fO) 

and by repeating this procedure n times one obtains 

V"|*) = |+) ® K"|^) + |-) ® XK"|^>. (f f ) 

Therefore 

h nl (12) 

= |+) ® e~ xt \yj) + |-) ® Xe~ tKt \^). 



It remains 'only' to identify X and to show that K is 
Hermitian. There is no general method allowing for find- 
ing solutions to quadratic equations like Eq. (6), yet 
judging from its structure it seems reasonable to try 
X = X)n£>il n "F By explicit substitution one in- 

fers that this is indeed a good guess as long as 

-S+^S 2 + \g\ 2 (n+l) >n 

Cn = , — — , n > 0. (13) 

g \n + 1 

The solution we have just found generalizes the Susskind- 
Glogower operator [8], in the sense that X — > (aat) _1 / 2 a t 
when S — > 0. As the result K = y/P + \g\ 2 (a^a + I B ) and 
it clearly is Hermitian. From Eq. (8) follows a = 
and in addition to that 

oo 

c(t) = e--*^e«e in " t (V'|n+l)HV) ! (14) 
n=0 

where Q n = ^ S 2 + \g\ 2 (n + 2) - s/P + \g\ 2 (n + 1). 

Note, we still have the freedom to choose the vector 
\ip). In particular, putting ~ |m), where \m) is a 
state with defined number of bosons, we have not only 
c(t) = 0, meaning p{t) is a steady state, but also c(t) = 
i.e. it exhibits its classical nature. 

In general, states alike the one from Eq. (5) are entan- 
gled, hence there is no one-to-one correspondence be- 
tween qubit-boson density operators and reduced 
density matrices Tred^) ( , F|) of the qubit [9], However, 
if \ip) is an eigenvector of X, i.e. X\ip) = \\ip), A e C, 
then |*) = (|+) + A|-)) ® \ip) is separable. 

Generalization. So far we have showed how to deter- 
mine initial preparations of the specific composite qubit- 
boson system which guarantee no energy flow between 
qubit and its environment (boson), i.e. the reduced qubit 
dynamics mimics pure dephasing evolution. We also in- 
vestigated conditions under which from this broad class 
of states one can chose the separable ones. 

Now, we are interested in finding answer to a more gen- 
eral question: How to prepare initial state(s) of the com- 
posite system (and determine their separability), which 
result in no 'information flow' between the subsystems, 
reflected in changes of a qubit observable A. Such states 
have very much in common with dephasing states and 
for that reason we will keep this terminology. 

Previously, we have conducted our analysis by mak- 
ing use of a very specific exactly solvable model. Cur- 
rently, we will show that neither solvability nor the form 
of the interaction in this model is a sine qua non con- 
dition for building 'dephasing states' for general qubit- 
environment models. 

Let A be a 2 x 2 Hermitian matrix-a given qubit observ- 
able. Our objective is to determine p = p(Q) such that 
(A(t)) = Tr(Ap(i)) remains constant during the evolu- 
tion. As before, by p{t) we denote the qubit reduced dy- 
namics, p(t) = Tr E [e _iH *|*)(*|e iH *], where |#) is the 
initial qubit-environment state and 

H = H Q ®lE + lQ®H E + Hi nt , (15) 
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with all the symbols having their usual meaning, stands 
for the total Hamiltonian. 

We begin with a very simple observation that in each 
moment of time t and for every complex 2x2 matrix 
A, the partial trace, Tr E (-), satisfies Tr [ATr E (£)(i))] = 
Tr[(A(g)I E )|*(i))(*(t)|]. By virtue of this relation 



(A(t)> = Tr[(A d ®I E )|Q(t))(Q(t)|], 
with \Q(t)) =e- ,Kt \n), where |fi) =U®I E |*), 



K = (U(8 I E )H(U g> I E 



it 



(16) 



(17) 



and U denotes the unitary matrix such that U'AU = Ad, 
with Ad = diag(A + , A_). The Kamiltonian K can always 
be written as 



K = |+)(+|®H + + |-)(-|®H. 
+ |+)(-|®V+|-)(+|®V t . 



(18) 



An explicit form of H± and V can easily be recovered 
when the Hamiltonians Hq, H e , and H; nt are provided. 
A priory we neither impose any physical restriction of 
their specification nor assume existence of symmetries in 
the total system (resulting e.g., in its solvability). 

Initially the composite system is assumed to be in the 
state g = where 



|*>=CV(|A + )®M + |A_>®X|V» 



(19) 



with |A±) = U + |±). \tp) is a freely chosen state of the 
environment. As before, one can redefine the state 
so that \tp) — > C^,\tp). We show that \Q) undergoes the 
following evolution 

|fi) -> e- lKt |fi) = |+) <g |^) + h) O X|&>, (20) 

provided that X solves the (operator) Riccati equation: 

XVX + XH + - H_X - V f = 0, (21) 

where \tpt) = e~' iK +*|^), for some operator K + . Riccati 
equations of the type (21) have been studied not only 
in the broad mathematical context [10] but also have re- 
cently been applied to investigations of open quantum 
systems [11]. Let us also notice here that Eq. (21) re- 
duces to the condition (6) when H is given by (1), as one 
may anticipate. 

Now, let K + = H + + VX, then from (21) immediately 
follows that XK + = H_X + V 1 , which leads to 

K|fi) = |+>®K + |V) + |-)«iXK + |V>. (22) 

As before one can easily justify the general formula: 

K"|f2) = |+)(g)K^) + |-)®XK^>, n>0, (23) 

from which the evolution generated by K is easily ob- 
tained to be 



\m) 



(24) 



Xe"' iK +*| 



Having (24) in place, one can formulate conditions un- 
der which (16) does not depend on time. A necessary 
condition for such evolution is for XV to be an observ- 
able. Because then K + is Hermitian and hence \ifit) evolve 
in unitary fashion. It is a matter of straightforward cal- 
culation to show that in this case Tr E (|^(0) has 
the dephasing structure (4), where 



a=(V#), c{t) = (^\e tK+t Xe 



-iK+t 



IV>>, 



(25) 



and therefore from (16) follows (A(i)) = aA + + (l — a)X_. 

Interestingly, the dephasing dynamics can also be de- 
signed by starting from any state orthogonal to (19). In- 
deed, such states are found to be of the form: 



i$) = a 



* (|A-> O |0> — |A+> <8) X+|0>) , (26) 



which can be verified by noticing that (^|$) = 0, for 
every \<p) € Ji^. In this case the dephasing dynamics 
takes the form 



o 



1 



c(t) = (0|e A 



£ XV 



provided K_ = H_ — V^X^ is Hermitian. For instance, in 
the Jaynes-Cummings mode l discussed in the preceding 
section we find that K_ = \/<5 2 + \g\ 2 a^a. 

Second example. As a second example we consider the 
fc-photon Rabi model [12], for which Hq = uxj z , H e = 
va) a read exactly as in the James-Cummings model (1) 
but the interaction is given by 



H int =<t x ® {g*a k +g{a' f ) k ). 



(28) 



This model not only generalizes the single mode case 
but also includes counter-rotating- wave terms [13], which 
make its analytical treatment much more complicated in 
compare with (1). Despite recent progress in finding it 
analytic solution [14], the problem remains highly non- 
trivial. 

Let us suppose that this time we want to find ini- 
tial state(s) of the composite system such that the x- 
componcnt of the qubit spin operator remains constant 
during the evolution i.e., S x = ^Tr (o~ x p(t)) = est. First, 
one needs to determine U which transforms o~ x to its di- 
agonal form. It is an easy task to do and the answer is 
U = (<r z + a x ) /a/2- In this case we have Ucr^U^ = a z , 
that is X± = ±1. Next, we transform H into K accord- 
ing to (17) and then recover H±, V from the decomposi- 
tion (18). Straightforward calculation shows 

H± = va ! a± [g*a k + g(a f ) k ) , V = wI E . (29) 
The corresponding Riccati equation (21) reads as follows 



uX 2 + XH, - H_X - wI E = 0. 



(30) 



We can solve this equation by introducing the generalize 
parity operator [15]: 



k oo 



(31) 



1=1 n=0 
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where \n,l) :— \kn + 1 — 1) and I < k. Such operator is 
both Hermitian and unitary, hence also Xj? = Ie- Note 
that for k = 1 this parity simplifies to the well know 
bosonic parity operator P = exp(z7ra^et). As one can 
verify X^a^X^ = — a k , thus X^H+Xk = H and therefore 
Xfe indeed solves (30). This fact may come as a surprise 
since (31) does not depend on any of the parameters is, 
ui, g. In this model the dephasing dynamics (25) and (27) 
are generated by K± = H± ± wX^ , respectively. 

Introducing projections P± = |(Ib ± Xj.), onto sub- 
spaces T-L± consist of states with defined parity (with 
respect to the generalized parity (31)) and taking into 
account both (5) and (26) we have 

l*e) = \ (1+) ® PeM + h) ® P-e|^» , e = ±1, (32) 

which are separable if € H±. 

Summary. The more conservation laws are present in 
the system (either classical or quantum) the better is 
our understanding of its behaviour and properties. In 
this Brief Report we have presented a method of finding 
time evolutions of quantum systems guaranteeing con- 
servation of certain quantities in the absence of related 
conservation laws. This is done by a proper choice of 
initial states. 

In order to built initial states resulting in the 



observable-conserving evolution one needs to solve the 
corresponding Riccati equation. Unfortunately, neither 
the form nor even the existence of such a solution can be 
taken for granted in general. Nevertheless, some useful 
criteria for solvability, applicable to the broad range of 
relevant physical systems, can be found in literature [16]. 

Preparation of the observable-conserving states or 
|$) would require highly sophisticated quantum engineer- 
ing. One could construct a state |x) = X\ip), then 'ten- 
sorize' the result with a qubit state |— ) ® \x), and finally 
superpose with the tensorized qubit-environment state 
|+) ® \ip). Clearly, the situation becomes simpler if the 
states are separable as it is the case for the broad class 
of problems when is an eigenstate of X. 

Despite obvious difficulties there is also at least one 
advantage of this procedure: It is performed only once 
at the beginning of the evolution-no dynamical control is 
needed to maintain the desired dynamics. We hope that 
together with continuous development of the quantum 
state engineering techniques the construction proposed 
in this work can become useful for applications. 
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